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Abstract 

Experimental data on the shape of hadronic momentum spectra are compared to theoretical pre- 
dictions in the context of calculations in the Modified Leading Log Approximation (MLLA), under the 
assumption of Local Parton Hadron Duality (LPHD) . 

Considered are experimental measurements at e + e~-colliders of the position of the maximum in 
the distribution of £ p = log(l/iEp), where x p — p/p D eam' The parameter £* is determined for various 
hadrons at various centre of mass energies. The dependence on the hadron type poses some interesting 
questions about the process of hadron-formation. The dependence of on the centre of mass energy 
is seen to be described adequately by perturbation theory. A quantitative check of LPHD + MLLA is 
possible by extracting a value of a s from an overall fit to the scaling behaviour of Q. 







1 Introduction 

During a few years of LEP running, a large amount 
of information has been collected on identified had- 
ron species in jets. The higher centre of mass energy 
of LEP, compared to past e + e~ accelerators, makes 
it easier to separate the behaviour of hadrons with 
a high momentum, which are correlated strongly 
to the primary quark, from those with a low mo- 
mentum, created mainly during fragmentation and 
hadronisation. 

Clear scaling violations have been observed in the 
shape of the charged particle distribution of x p — 
p/E-fo ea;nl as a function of the centre of mass energy. 
The strong coupling constant has been extracted 
from these scaling violations, using the behaviour at 
high momenta: 0.2 < x p < 0.7 JlJ]. Due to the larger 
statistical errors this is not possible for individually 
identified hadron species. 

However, low momentum data for specific types 
of hadrons may be used to study the properties of 
jet-evolution and hadron-formation in the context 
of the LPHD hypothesis and MLLA calculations of 
parton spectra ^MQ,|j. The assumption of 'Local 
Parton Hadron Duality' (LPHD) states that a cal- 
culated spectrum for 'partons' in a 'parton shower' 
can be related to the spectrum of real hadrons by 
simple normalisation constants. These constants 
have to be determined by experiment. A second 
assumption is that the low momentum part of the 
spectrum is not influenced in a significant way by 
hadrons that are correlated to the primary quark. 

Calculations of the parton spectra in the 'Modi- 
fied Leading Log Approximation' (MLLA) take into 
account next-to-leading logarithms in a consistent 
fashion. The physical mechanism relevant to these 
next-to-leading logarithms is the coherent emission 
of soft gluons inside a jet, leading to an angular or- 
dering and an effective transverse momentum cutoff 
for the partons. Parton jets develop through re- 
peated parton splittings, resulting in an increase of 
the multiplicity at lower momenta. The interplay 
of coherent emission of gluons and the creation of 
hadrons causes this spectrum to be cut off at very 
low momenta. Calculations predict the shape of 
the distribution of £ p = log(l/a;p). The resulting 
'hump-backed' distribution is nearly gaussian. As 
an example, the Monte Carlo spectrum of the A 
baryon in Z° decays can be seen in figures [l] and 
||. In the following, the maximum £* in the £ p 
distribution will be determined for various types 
of hadrons and at various centre of mass energies. 
Subsequently a comparison is made with the theo- 
retical calculations. 
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Figure 1: The dependence of the A spectrum on 
the species of its 'parent' particle in Z° decays simu- 
lated by JETSET. Particles produced directly in the 
fragmentation process have the Lund-model string 
as 'parent'. The scales of both x p and £ p are given 
for comparison. 
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Figure 2: The dependence of the A spectrum on 
the flavour of the primary quarks in Z° decays at 
LEP, as predicted by JETSET. Note that there is not 
a clear separation into a region of high momentum 
(correlated to primary quarks) and another of low 
momentum (related to jet fragmentation). 
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charged 


7T± 


77° 


K± 


K° 


P 


A 




ARGUS 
CLEO 


9.98 
10.49 




2.32 ±0.022 
2.30 ± 0.08 


2.42 ± 0.08 


1.64 ±0.05 


1.72 ±0.04 
1.68 ±0.17 


1.63 ±0.07 
1.67 ± 0.10 


1.44 ± 0.034 
1.47 ± 0.08 


1.32 ± 0.11 


TASSO 
JADE 


14 
14 


2.45 ±0.05 


2.66 ±0.06 


2.68 ± 0.27 


1.91 ± 0.13 


1.65 ±0.13 


1.80 ± 0.15 






TASSO 
JADE 


22 
22.5 


2.74 ±0.06 


2.99 ±0.09 


2.84 ± 0.30 


2.41 ± 0.23 


2.47 ±0.8 


2.14 ± 0.27 


1.75 ±0.50 




HRS 

TPC/27 


29 
29 




> 3.3 
3.00 ± 0.05 




2.28 ±0.42 
2.14 ± 0.05 


2.25 ±0.30 
1.98 ±0.09 


2.13 ±0.11 


1.96 ±0.12 




TASSO 
JADE 


34 
35 


3.00 ±0.05 


3.17 ±0.05 


3.24 ± 0.05 


2.35 ±0.20 


2.32 ±0.10 


2.24 ± 0.09 


2.15 ±0.15 




TASSO 
JADE 


44 
44 




3.37 ± 0.07 


3.52 ± 0.08 


<3.5 




<3.5 






TOPAZ 


58 


3.42 ± 0.04 


3.51 ± 0.07 




2.83 ±0.10 




2.64 ± 0.06 






OPAL 

DELPHI 

L3 


91.2 
91.2 
91.2 


3.60 ± 0.04 
3.71 ± 0.05 


3.79 ± 0.04 


3.96 ± 0.13 


2.73 ± 0.06 


2.87 ±0.07 
2.62 ±0.11 
2.89 ±0.05 


2.95 ± 0.09 


2.78 ±0.08 
2.81 ± 0.04 


2.55 ± 0.20 



Table 1: The values of £* for various hadrons as determined from the momentum spectra as measured 
by e + e _ experiments at different values of i/s. For the rj meson, the L3 experiment reported a value 
C = 2.60 ±0.15. 



Figure |] shows that LPHD is not an obvious as- 
sumption: in the JETSET Monte Carlo || , the spec- 
trum of e.g. a A baryon depends on its 'parent' par- 
ticle. The LPHD assumption states that the sum of 
these spectra is proportional to the spectrum as cal- 
culated for a parton shower with a correctly chosen 
energy cut-off. 

Moreover, the JETSET momentum spectrum of 
the A baryon (figure ||) shows a strong dependence 
on the flavour of the primary quarks. For the heav- 
ier quarks, s, c and b, the flavour dependence does 
not only manifest itself at high momenta, but also 
affects the distribution at low momenta. This will 
cause shifts of £* that depend on the jet-flavour. 

The spectrum of all charged particles at LEP has 
been fitted to the MLLA distribution, with free nor- 
malisation factors for pions, kaons and protons, de- 
pending on the centre of mass energy |?],|| ■ It was 
quite surprising Q that the MLLA functions can 
also fit the high momentum part in the data for pi- 
ons. This is probably a coincidence, since the heav- 
ier K° meson has a spectrum that can not be fitted 
as nicely at large x p . 

One way to see whether the observed shape of the 
spectrum is really due to coherent gluon emission 
is to compare experimental data to predictions of 
the JETSET Monte Carlo program, with the coher- 
ence (angular ordering) either turned on or turned 
off. The large dependence on the primary quark 
flavours and the fact that coherence in JETSET is not 
really necessary to fit the experimental data has led 
some to the conclusion H that coherence can not 



be demonstrated on the basis of hadron momentum 
spectra. 

Here an attempt will be made to infer more in- 
formation using the properties of identified hadrons, 
concentrating on the observable £* . This is done for 
data from LEP and from various other e + e~ collid- 
ers. Subsequently, the available results for Q_aie 
compared to the LPHD + MLLA approach 



Experimental values of £! 
identified hadrons 



for 



At LEP, analyses of explicitly identified hadrons 
have been performe d for the vt* [[To) , the n° , 

the i] pljilf the pro- 



the Kg 



the K ' JU^ 

ton Jh} and the A [Lj.|lJ|Hi[[ljJ . At lower statistics 
data is also available for the S~ and Q~ baryons, 
as well as the S* and S* 

The values of £* correspond to low momenta, 
where the dependence on the primary quark fla- 
vours is expected to be small. To first approx- 
imation the distribution in £ p is gaussian, but 
a distorted gaussian fits the calculated spectrum 
more accurately [|l8]|| ■ The statistics of identified 
hadrons is not sufficient to extract the width and 
the two distortion parameters for more than a few 
experiments. 

However, since these other parameters depend 
more strongly on the high momentum tail of the 
distribution, they can be expected to depend more 
strongly on the event flavour, which could bias 
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Figure 3: The maxima £* of the distributions in £ = — logx p for different hadrons as measured at LEP 
(table |) are shown in figure (a). Two different determinations of the corresponding scale Qo are shown 
in figures (b) and (c). The calculated dependence for A e g? = 50 MeV leads to figure (b), where one would 
naively expect Qo to be proportional to the hadron mass. Figure (c) gives the values of Qo as determined 
from the fit of the dependence of £* on E^ eam . A (large) common error on the different values of Qo has 
been neglected in both figures (b) and (c). This causes the vertical scales to be rather arbitrary, but the 
relative position of the points remains significant. 



the results. We will suppose that this violation of 
LPHD can be avoided adequately by concentrating 
on £ * , a property of the low-momentum part of the 
spectrum. 

In this paper £* is simply defined as the maximum 
of the distribution. The values were determined by 
fitting the £ p distribution in a limited range (~ ±1) 
of £ p around the maximum, to a gaussian distribu- 
tion. A systematic error was estimated by changing 
the fitted range and by fitting the same range to a 
polynomial function. The statistical and systematic 
errors were added in quadrature. The systematic 
error often dominated. 

For the most accurately determined experimen- 
tal spectra, those of identified charged hadrons at 
OPAL jnj, a more complicated function is neces- 
sary. In this case a polynomial of order 4 or 5 was 
used to fit either the partial cross section, or its log- 
arithmQ These polynomials were defined such that 
the maximum was one of the free parameters in the 
fit. In this way the MINOS procedure of MINUIT 
p7[ can automatically determine an interpolation 
error due to variations of the other free parame- 
ters of the polynomial. The values found for the 
maxima in OPAL's spectra are different from those 
reported by OPAL, and the estimated errors are 
slightly larger. 



The results are given in table 0. Use was made of 



spectra published by CLEO 
21), HRS 



|2 



9|, ARGUS |20[,TPC 
JAD E p| TOPAZ 
§[l||l5]g and L3 



J§, TASSO [| 
DELPHI HQ, OPAL 
1 8| Jl4] , |TT[ | . Most of the recent publications reported 
values of £*, but they did not all use exactly the 
same definition of it. While the numbers given by 
L3 were compatible with the method used here, 
OPAL and DELPHI used slightly different defini- 
tions and for consistency the values of £* are those 
determined by the author from the published ex- 
perimental spectra. The spectra from TOPAZ were 
not available in |p5| , only values of £*; likewise for 
the recent paper by L3 Q. For most experiments 
at lower centre of mass energies the authors did not 
report the distribution of £ p or x v , but used the en- 
ergy to define xe = 2-E/-E Deam , or they divided by 
(3 = v/c to obtain the 'scaling cross section'. Many 
different normalisations have been used, but fortu- 
nately this is unimportant for a determination £*. 
From the experiments (ARGUS and CLEO) near 
the T resonances results are given for the 'contin- 
uum', and for the resonances itself, where the con- 
tributions from the continuum are subtracted. The 
values in table |l| refer to data from the continuum 
above the T resonance. In the data from ARGUS 
for charged pions and protons, the contributions 
from decays of Kg and A were subtracted. 



1 Fitting the logarithm to a 4th order polynomial is equiv- 
alent to using a distorted gaussian. 
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and A baryons. The curves for kaons and protons 
happen to lie on top of each other. 



3 Predictions for £ 



The LPHD + MLLA calculations J| of the dis- 
tribution of £ p depend on an effective QCD scale 
^QCD ~ ^eff an< ^ 011 a transverse momentum cut- 
off Q hi the evolution of the parton cascade. The 
value of £* is calculated to be a nearly linear func- 
tion of logQo, log A eff and log#b eam . 

The calculated dependence on the centre of mass 
energy is as follows [|: 



Y 



- + tJC/Y - C/Y + 0(Y- 3 / 2 ) 
F(X), 



where F(0) = and f\ 
Y = log(£ beam /A eff ), 



A = log(Q /A eff ), (2) 



and A e Q is an effective QCD scale, while only the 
cutoff scale Qq depends on the hadron type. The 
constant C is calculated to be 



where 



HJV c /3 + 2n//(3JV*), 
lL/V c /3-2n//3. 



(3) 



(4) 
(5) 



Figure 4: The experimental values of £* versus ^fs 
for various identified hadrons. The curves are from 
the combined fit of A e g: and Ffhadron) to the pre- 
dicted dependence (equations (0) and (|TT|)). From 
top to bottom, the curves correspond to pions (n^ , 
7T°), all charged particles, kaons (K^, K°), protons 1S: 



N c = 3 is the number of colours and rif = 3 is the 
active number of quark flavours in the fragmenta- 
tion process. For n/ = 3 one finds C = 0.2915. 

It is important to note that F does not depend 
on Y, and that the first part of equation ([!]) is in- 
dependent of Qq. This predicted behaviour can be 
checked by comparing spectra of different identified 
particles. In the available momentum range this 
leads to a nearly linear dependence of £* on Y. 

In ref. || various graphs show the results of nu- 
merical calculations for the dependence of various 
parameters on E^ eam and Qq ^ ^eff' Fig ure 4 
of ref. H shows the dependence of £* on Qq for 
Agg = 150 MeV, and at a number of beam en- 
ergies. From this it is possible to extract F after 
subtracting the Y dependent part of equation (|l|). 
For an interpretation of the experimental values of 
£* it was useful to fit both this function F(X), and 
its inverse X(F) to polynomials. The result for F(X) 



F(X) — —1.46 • A + 0.207 • A 2 ± 0.06 (6) 

while the inverse of F is described by: 

X(F) = —0.614 • F + 0.153 • F 2 ± 0.06 (7) 

The given errors denote the maximum deviation 
from the distribution plotted in ref. ||), in the 
ranges -2 < F < and < A < 2. 

4 The meaning of A e g> 

The QCD scale A e g? [^6| is related to the (running) 
strong coupling constant by: 



2tt 



1 

bY 



^og^beamMcff)' 



(8) 



(1) 



2 This detail 
where equation 



is not completely clear in reference 
is given as a function of Y — X, but in the 
approximation that A = 0. 
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Figure 5: The values of Q' from table |] versus the 
mass of the hadron. 



This makes it possible to express the scale depen- 
dence of eq. (0) in terms of a s : 



BY 
dY 2 ' 



1 

8° 



-ab 



2irN c 

3/2 



o 



16VK 



■Ola 



7/2 



(9) 
(10) 



The y/al term is the next-to-leading (MLLA) cor- 
rection to the leading term. The next-to-next-to- 
leading term 0(a s ) vanishes |p6j . 

Although this sounds impressive, it is still not 
enough for an unambiguous definition of A e g or, 
equivalently, of a s (E) at a well-defined energy scale 
E. An uncertainty of e.g. a factor 2 in the energy 
scale leads to a s (2E) instead of a s (E). This induces 
a correction in equation (0) of 0{aV 2 {E)) 1 which 
shows that for a better definition of the energy scale 
it is necessary to know the 0{a Z J 2 {E)) correction. 
This is the reason for talking about the scale A e g, 
in stead of a more well-defined scale like A^r. 

If we now neglect the dependence of a s {E) on E, 
£* becomes a linear function of log(£ , ^ eam ), where 
the slope is given by equation (||). 



5 Dependence of £* on ■ s fs 

A more quantitative comparison is possible for the 
scaling behaviour of £* , as parametrised by the ef- 
fective QCD scale A e g-. The dependence on y^s = 
2_Eu eam of the experimental values of £* from ta- 
ble ]l| is presented in figure ||. 



As a start, the experimental values can be fit- 
ted to separate first- or second order polynomials 
in ^>gE^ eaTa for each type of particle. The qual- 
ity of these fits provides an internal check of the 
estimated experimental errors. Table |^ gives the 
slopes of fitted straight lines and the values of x 2 
for first and second order polynomials. The quality 
of the fits is acceptable, but the values of x 2 when 
fitting so many free parameters could indicate that 
the experimental errors have been underestimated 
somewhat. It seems that the measured slopes are 
not quite equal to one other, as was predicted by 
LPHD + MLLA. This is due to some dependence 
on the flavour of the particles and of the primary 
quarks, violating the LPHD assumption slightly. 

For a comparison with other publications it is 
also interesting to fit the data for each hadron 
type to equation ([!]), under the assumption that 
Q' = A cff = Q , so that F(X) = F(0) = 0. The 
energy scale Q' is different from Qq as determined 
in the previous section, and it is also different from 
Agg. In many past publications the notation l A e g 
is also used for Q' Q . This is confusing, since Q' is not 
expected to be the same for different hadrons, while 
the QCD scale A e Q should be independent^]. The 
quality of this fit with A c g- = Qq = Q' is slightly 
worse, with an overall \ 2 — 2.41. The largest con- 
tributions to x 2 come from the kaons, that also have 
a different slope bh from the other hadrons. Look- 
ing carefully at figure ^ one might also suspect the 
measurements done at the Z° and T. 

The fitted values of Q' are not equal to the corre- 
sponding values of Qo shown in figure [?] The values 
of Q' are quite near to the particle mass for mesons, 
but for the baryons Q' is significantly smaller. Fig- 
ure U shows the dependence of Q' Q on the hadron 
mass. No clear linear dependence is observed. 

To check how well the data is described by the 
calculations, a simultaneous fit of equation ([!]) was 
done to all the data from table [l], where A q q is a 
universal scale, while F is allowed to have differ- 
ent values for hadrons, pions, kaons and the two 
baryons. 

Such a six parameter fit was performed by min- 
imising x 2 using the MINUIT |27| program, with 
zero as starting values for F^. There is a strong 
correlation in the fit of log A„q and F, the average 
of Fh- Therefore it is practical to fit F and differ- 
ences Fh — F. After this, the fit converges properly, 



3 In principle one might expect that different hadrons 
from different decay chains could introduce a small depen- 
dence of an effective value of nj on the type of hadron. How- 
ever, equation (^|) shows that variations of n j from 3 to 5 can 
only change the slope by a few percent. This is well within 
the other uncertainties. 
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particle 


ah 


b h 


X 2 /n.d.i. 


Qo 


xVn.d.f. 


points 


charged 


0.76 ±0.12 


0.641 ± 0.030 


1.36 


0.212 ±0.007 


1.10 


6 


7T° 

7T° 


0.79 ± 0.06 
0.78 ±0.20 
0.78 ± 0.05 


0.667 ±0.018 
0.705 ± 0.059 
0.679 ±0.016 


1.23 
0.42 
1.49 


0.179 ±0.005 
0.147 ±0.008 
0.174 ±0.004 


1.53 
0.70 
1.91 


9 
6 
15 


K± 
K° 


0.47 ±0.11 
0.49 ±0.09 
0.48 ±0.06 


0.512 ±0.032 
0.517 ±0.025 
0.516 ±0.019 


2.17 
2.12 
1.93 


0.645 ± 0.027 
0.620 ± 0.023 
0.631 ±0.018 


4.48 
5.49 
4.78 


8 
10 
18 


P 
A 

p, A 


0.27 ±0.12 
0.02 ±0.07 
0.09 ± 0.07 


0.580 ± 0.036 
0.613 ± 0.021 
0.609 ±0.019 


0.37 
0.30 
1.29 


0.625 ± 0.030 
0.759 ± 0.025 
0.712 ±0.020 


0.93 
1.07 
1.69 


8 
7 
15 


total x 2 /n.d.f. 


1.21 


2.41 


54 



Table 2: Result of a fit of the experimental values of £* in table [T] to a linear function £* — ah + bh ■ 
\og(E cms /GeV) for various identified hadrons. The scale Q' is the result of a second fit to the function 
= \Y + yf(JY - C, where Y = log(.E beam /Q ) and C = 0.2915. The fits improve slightly when data 
at the T are excluded. To compare the overall quality of the fits to the two functions, a total x 2 /n.d.f. 
is given, based on the fits for n , 7r°, K^, K°, p, and A. 



but the correlation coefficient of logA e g- and F is 
still 100%. 

The strong correlations can be understood from 
equation (0) and figure ||: a change in log A trans- 
lates the curves in the direction of the x-axis, while 
a change of F causes a translation along the j/-axis. 
The fact that the curves are so straight causes the 
strong correlations. The reason that the fit still con- 
verges is that the slope of the curves varies slowly 
with changing A c g\ 

To see if this causes numerical instability or even 
inaccuracy, it was subsequently checked that the 
correlation in the fit can be reduced to a more 
acceptable 74% by a change of variable to F' — 
F — r ■ (log A — log Ao), where r = 0.62 is the aver- 
age value of dt;*/dY as given by equation ([!]) and 
the constant Ao ~ A e g-. The fit in this parametri- 
sation converges to the same values of x 2 and A e g 
as when fitting F itself. 

The fit has 47 degrees of freedom and leads to 
X 2 /47 = 2.15. This value of x 2 should be compared 
to the values of the 'total' x 2 / n -d.f in table ||, that 
are also larger than 1. 

The resulting parameters and their statistical ac- 
curacy (after applying the scale factor \J x 2 /n.d.f. 
1.17 to all the experimental errors) are: 



A cff 

=4> a s (m z o) 
F 

F^-F 

^charged — * 
F K -F 
F p -F 
Fa 



0.052 GeV 
0.103 ±0.022 

+0.75 
-1.1 

±0.021 



F 



-1.29 
0.54 
0.42 
-0.28 
-0.28 
-0.41 



±0.027 
±0.027 
±0.039 
± 0.030 



(11) 



The errors here are purely statistical, as given by 
the MINOS method in MINUIT @, which takes 
into account non-linearities and correlations be- 
tween parameters. The correlation coefficient of 
A e g and F is 100%, while the values of Fh — F 
are constrained to have a zero sum. 

Figure ^ shows that the fitted functions are very 
near to straight lines and that they follow the data 
points rather well. Addition of a free 0(Y 3 ^ 2 ) term 
as in equation ([!]) does not improve the fit signifi- 
cantly. 

So what should one conclude from the high val- 
ues of x 2 ? It seems that the general behaviour is 
described well by LPHD ± MLLA, but that equa- 
tion (Q) differs significantly from the data for each 
hadron. This is probably a signal of the breakdown 
of LPHD at this level of accuracy. 

An alternative way of fitting the data is based on 
equation ([)]) with a fixed value of a s . Allowing a s 
to 'run' somewhat leads to the following parametri- 
sation: 



= Poh + Pi (log ^beam/ GeV _ lo S i5 ) 
±ip 2 (log£ beam /GeV-logl5) 2 , 



(12) 



where the quadratic term is centred at 15 GeV. The 
constant term po/i is different for each hadron type. 
Equations (0) and ( |l0| ) relate the values of pi and 
P2 to the strong coupling constant a s . This fit has 
no problems with correlated parameters and it con- 
verges to: 



Pi 

P2 



0.613 ±0.015 
-0.012 ±0.026 



+0.035 
-0.030 



0.122 
0.14 ±0.16 



(13) 



6 



where a s is given at a scale of 15 GeV. Again, 
the given errors are purely statistical. The value 
of A eff from equation (jl^) gives a s = 0.123 io!o3i 
at the scale 15 GeV, which is equal to the value 
derived from p\. Within its large error, pi is also 
consistent. The values of x 2 / n -d.f are comparable 
to those in the previous fit. The differences Fh — F 
can be determined from poh arid they are consistent 
with those given in eq. (pT|). 



6 Dependence of 
hadron mass 



on the 



The dependence of £* on the hadron mass as deter- 
mined from LEP data is shown in figure ||(a). As 
has been said, this dependence is not predicted by 
the LPHD + MLLA approach, but it is interesting 
to look at the relation between the scale Qq and the 
mass and flavour of the hadron. Using equations (|]) 
and (0) and the determined value of A p p- = 50 MeV, 
one can convert the values of £* (eq. [11 ) to values of 
Qo/Agfj. The result of this is shown in figure ||(b). 
The error bars in this figure do not incorporate the 
uncertainty of A e g-. This means that the absolute 
scale should not be taken too seriously. 

Naively one could expect that the cutoff scale 
Qo grows proportionally to the mass of the hadron. 
The most naive guess would be a linear dependence 
of the form Q w A cff + m hadron . Figure |(b) 
shows that this is not correct and that at least a 
separate treatment of mesons and baryons is neces- 
sary. 

Figure ||(c) also shows Qo, but now determined 
from the fitted Fh in equation ( p"i"| ) . Again the errors 
given do not include the very large overall error on 
F and they are only significant when comparing the 
different points relative to each other. Taking this 
into account, the result is comparable to that of 
figure |(b). 

7 Summary 

The data of LEP are giving a wealth of informa- 
tion about the fragmentation of quarks and gluons, 
as well as the mechanisms at play in the hadroni- 
sation. Although Monte Carlo models like JETSET 
and HERWIG can describe these data very accurately, 
it is a valid question whether these models with all 
their free parameters lead to a better understand- 
ing of the physics behind the formation of hadrons 
in a jet, or are merely to a better parametrisation. 

The 'LPHD + MLLA' approach §§J|[| is an in- 
teresting attempt to understand some properties of 



multihadron production in jets in terms of pertur- 
bative QCD. This approach tries to take perturba- 
tive QCD as near as possible to the limits posed 
by confinement. It gives predictions for properties 
of momentum spectra that can be compared with 
experimental data, but concentrates on observables 
that can be described without too many unknown 
free parameters. 

The experimental values of the parameter £* for 
various hadrons and at various centre of mass ener- 
gies was determined from the published momentum 
spectra. Subsequently, an investigation was made 
of the dependence of £* on both the centre of mass 
energy and the mass of the identified hadrons, and 
this could be compared with the theoretical predic- 
tions. 

The dependence of £* on the mass and flavour of 
the identified hadron raises some interesting ques- 
tions that have not yet been answered in the context 
of LPHD + MLLA and the description of fragmen- 
tation by truncated parton cascades. 

An attempt was made to clear up some of the con- 
fusion in the definition of A e g. The author would 
like to stress that it is a common but confusing prac- 
tice to use the symbol A„g for the other scale, here 
called Q' . This is because the Q' is not expected to 
be independent of the type of hadron, and should 
therefore not be interpreted as Aqqjj). 

The price paid for a more consistent definition of 
A e g is the explicit introduction of the constants F, 
that are different for different hadrons. The advan- 
tage is that A e g can now be interpreted as a QCD 
scale Aqqjj), and that it is then possible to make a 
more quantitative comparison with theory. 

The dependence of £* on the centre of mass en- 
ergy is described adequately by the MLLA calcu- 
lations. A value of a s has now been extracted 
from these data, and its value is consistent with 
the many accurate measurements. It is perhaps 
striking that data at very low momenta (at LEP, 
p » 0.83(2. 27)GeV when £ p w 3(4)), that are near 
to the region of phasespace where confinement oc- 
curs, can not only be described qualitatively by per- 
turbative QCD, but can even be used to extract a 
consistent value of the strong coupling constant. 

It is a success of LPHD + MLLA that the ex- 
tracted value of a s is correct. However, it is clear 
that this is not a good way to determine a s ac- 
curately: it is nearly impossible to make a proper 
estimate of the systematic errors. The quality of 
the overall fit suggests that the picture of LPHD is 
starting to break down. 
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